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'. Abstract 



^ ■ We give a comprehensive review of the renormahzation group method for global 

and asymptotic analysis, putting an emphasis on the relevance to the classical theory 
\ of envelopes and the existence of invariant manifolds of the dynamics under consid- 

eration. We clarify that an essential point of the method is to convert the problem 
from solving differential equations to obtaining suitable initial (or boundary) condi- 
tions. We mention that the notion of envelopes is also useful for constructing global 
• and asymptotic behavior of wave functions of quantum systems such as the ones 

. with the quartic potential or double-well potential. 

o> 
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1 Introduction 



, The renormalization group (RG) equations have a peculiar power to improve the global 
nature of functions obtained in the perturbation theory in quantum field theory The 
RG equations represent the fact that a physical quantity 0{p,a, /i) should not depend on 
^ ' the renormalization point /i, 



dOjp, a; /i) ^ Q 

Such a floating renormalization point was first introduced by Gell-Mann and Low in the 
celebrated paper |^. We remark that the renormalization point fi plays a role of the initial 
point and the renormalization condition for physical quantities such as coupling constants 
at the energy scale fi may be viewed as setting initial values of these quantities. 

Recently, the present author has indicated that the RG equation a la Gell-Mann-Low 
can be identified as the envelope equation]^;, it was shown that the notion of envelopes is 
useful for improving the perturbative expansions appearing in quantum field theory. This 
identification was realized through an examination of the RG method of Chen, Goldenfeld 
and Oono for global analysis of differential equations 0; they applied the RG equation to 

^Talk presented at RIMS (Kyoto) Workshop on Geometrical Methods for Asymptotic Analysis 
1997.5.20 - 5.23 
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numerous problems and found that the RG equation gives slow dynamics of the system 
in question. The new point of their method is to utilize secular terms which usually 
appear when perturbation theory is applied to differential equations; this is in contrast 
to all previous methods, which are formulated on the principle to avoid the appearance 
of secular terms. 



Their method was reformulated on the basis of the classical theory of envelopes[^, ^. 
It was demonstrated that owing to the very envelope equation, the functions constructed 
from the solutions in the perturbation theory certainly satisfies the differential equation in 
question uniformly up to the order with which local solutions is constructed. It was also 
shown 1^ in a most general setting that the RG method gives a reduction of dynamics 
and clarified that there is a correspondence between the RG method and the reductive 
perturbation method]^ : Some interesting examples were also worked out in this method, 
such as the forced Duffing, the Lotka-Volterra and the Lorenz equations; the first example 
showed that the method is applicable to not only automonous but also non-automonous 
equations, the second one showed that the method gives phase equations and the last 
one showed that center manifolds of the dynamics can be extracted in the method. It is 
noteworthy that the notion of functional self- similarity (FSS) extracted as the essence of 
the RG in is only applicable to autonomous equations; see Q for an application of the 
notion of FSS for deducing phase equations. The above fact suggests that the notion of 
envelopes better represents the underlying mathematics of the powerfulness of the RG. 
It was recognized also that when the unperturbed equation has neutrally stable solutions, 
the RG method works well. It implies that the method is applicable when the dynamical 
system under consideration has invariant manifolds and useful to extract the manifolds 
and the dynamics on the manifolds. It was shown that global and asymptotic behavior 
is obtained even for discrete systems by constructing "envelopes" when the system has 
neutrally stable solution of the unperturbed equations [0. 



Subsequently the classical theory of envelopes was applied for getting asymptotic be- 
haviors of wave functions in quantum mechanics\10\. This is an optimized perturbation 
theory in which the perturbation theory is combined with a variational method. The key 
ingredient is to construct an envelope of a set of perturbative wave functions. This leads 
to a condition similar to that obtained from the principle of minimal sensitivity[|ll|]. Ap- 
plications of the method to quantum anharmonic oscillator and the double well potential 
show that uniformly valid wave functions with correct asymptotic behavior are obtained 
in the first-order optimized perturbation even for strong couplings. 



In this report, we will give only basic ingredients of the RG method based on our 
formulation putting an emphasis to the relevance to envelopes; for detailed account of the 
method and various applications, please refer to Ref.'sQ, H, ^, |, |10|- In this report, we 
also clarify that the RG method is most lucidly formulated by noting that the method 
converts the problem from solving differential equations to obtaining suitable initial (or 
boundary) conditions as does the usual RG. In this report, we shall not discuss on partial 
differential equations due to the lack of space. Please refer to as well as p, ^, |^ for 
this subject. 
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2 A short review of the classical theory of envelopes 



We here give a brief review of the theory of envelopes. Although the theory can be 
formulated in higher dimensions we consider here envelope curves, for simplicity. 

Let {Cr}r be a family of curves parametrized by r in the x-y plane; here C^- is rep- 
resented by the equation F{x,y,T) = 0. We suppose that {Ct-}^ has the envelope E, 
which is represented by the equation G{x,y) = 0. The problem is to obtain G{x,y) from 
F{x,y,T). 

Now let E and a curve Cj-o have the common tangent line at {x,y) = {xo,yo), i.e., 
(xo, yo) is the point of tangency. Then Xq and yo are functions of Tq; Xq = 0(ro), yo = '?/'(to), 
and of course G{xo,yo) = 0. Conversely, for each point (xci/o) on E, there exists a 
parameter tq. So we can reduce the problem to get tq as a function of {xo,yo); then 
G{x,y) is obtained as F{x,y,T{x,y)) = G{x,y)^ T"o(a^05?/o) can be obtained as follows. 

Since the tangent line of E at (xo,?/o) is perpendicular to the normal direction of 
F{x, y,T) = at the same point, one has 

Fxixo, yo, T-o)0'(ro) + Fy{xo, yo, ro)^'(ro) = 0. (2.1) 

On the other hand, differentiating F(x(ro), y(ro), tq) = with respect to tq, one also has 

Fxixo,yo,To)(l)'iTo) + Fy{xo,yo,ro)ip'{To) + F^^{xo,yo,To) = 0. (2.2) 

Combining the last two equations, we have 

Fro [Xo, yo, To) = ^ = 0. (2.3) 

OTo 

This is the basic equation of the theory of envelopes; we call this type of equation envelope 
equations and also RG/E equation where RG and E stand for renormalization group and 
envelope, respectively, becasue the RG equation Eq(l.l) has the same form as Eq.(2.3). 
One can thus eliminate the parameter tq to get a relation between xo and yo, 

G{x,y)=F{x,y,To{x,y)) = Q, (2.4) 

with the replacement {xo,yo) {x,y). G{x,y) is called the discriminant of F{x,y,t). 

Comments are in order here: (i) When the family of curves is given by the function 
y = f{x,T), the condition Eq.(2.3) is reduced to df/dro = 0; the envelope is given 
by 2/ = f{x,To{x)). (ii) The equation G{x,y) = may give not only the envelope E 
but also a set of singularities of the curves {C^jr- This is because the condition that 
dF/dx = dF/dy = is also compatible with Eq. (2.3). 

^ Since there is a relation G{xo, yo) = between xq and yo, to is actually a function of xo or yo- 
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3 The RG method; a simplest example 



In this section, using a simplest example we show how the RG method works for obtaining 
global and asymptotic behavior of solutions of differential equations. We shall present the 
method so that the reader will readily see that the notion of envelopes is intrinsically 
related to the method. We shall emphasize that an essential point of the method is 
tuning the initial condition at an arbitrary time to perturbatively along with solving the 
perturbative equations successively. One will see that the reasoning for various steps in 
the procedure and the underlying picture are quite different from the original ones given 
in 0. We believe, however, that the present formulation emphasizing the role of initial 
conditions and the relevance to envelopes of perturbative local solutions straightens the 
original argument, and is the most comprehensive one.Q 



Let us take the following simplest example to show our method: 

d'^x dx 



where e is supposed to be small. The solution to Eq.(3.1) reads 



xit) = Aexp(--t) sin(Wl - jt + e), (3.2) 



where A and 6 are constants. 



Now, let us obtain the solution around the initial time t = to in a perturbative way, 
expanding 

x(t,to) = xo(t,to) + exi{t,to) + e^X2{t,to) + (3.3) 
where Xn(t,tQ) {n = 0, 1,2...) satisfy 

Xo+ Xo = 0, Xn+l + Xn+l = "Xn- (3.4) 

The initial condition may be specified by 

x{to,to) = W{to). (3.5) 

We suppose that the initial value W(to) is always on an exact solution of Eq.(3.1) for any 
to- We also expand the initial value W{to); 

W{to) = Woito) + eWM + e^W^ito) + (3.6) 

and the terms VFj(to) will be determined so that the perturbative solutions around different 
initial times to have an envelope. Hence the initial value W{t) thus constructed will give 
the (approximate but) global solution of the equation. 

■^Tlic following is even a refinment of the argument given in ||, |^ where the fact that the RG method 
is a theory manipulating initial conditions were not fully recognized. 
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Let us perform the above program. The lowest solution may be given by 

xo{t, to) = A{to) sm{t + 9{to)), (3.7) 

where we have made it explicit that the constants A and 9 may depend on the initial time 
to. The initial value W{to) as a function of to is specified as 

Wo{to) = xo{to, to) = A{to) sin(to + 9{to)). (3.8) 

We remark that Eq.(3.7) is a neutrally stable solution; with the perturbation e 7^ the 
constants A and 9 may move slowly. We shall see that the envelope equation gives the 
equations describing the slow motion of A and 9. 

The first order equation now reads xi + xi — —Acos{t + 9), and we choose the 
solution in the following form, 

xi{t, to)^~-{t- to) sm{t + 9), (3.9) 

which means that the first order initial value Wi{to) = so that the lowest order value 
Wo{to) approximates the exact value as closely as possible. Similarly, the second order 
solution may be given by 

Mt) ^^{{t- to)^ sin(t + 9)-{t- to) cos(t + 9)}, (3.10) 
8 

thus W2{to) = again for the present linear equation. 

It should be noted that the secular terms have appeared in the higher order terms, 
which arc absent in the exact solution and invalidates the perturbation theory for t far 
from to- However, with the very existence of the secular terms, we could make Wi{to) 
{i = 1, 2) vanish and W{to) = Wo{to) up to the third order. 

Collecting the terms, we have 

x{t,to) = Asm{t + 9)-e—{t-to)sm{t + 9) 

+^^^{{t - to)^sm{t + 9)-{t- to) cos(t + 9)}, (3.11) 
8 

and more importantly 

W{to) = Woito) = A{to) sin(to + 9{to)) (3.12) 

up to O(e^). We remark that W{to) describing the solution is parametrized by possibly 
slowly moving variable A{to) and 4>{to) =to + 9 (to) in a definite way. 

Now we have a family of curves {Cto}t(, given by functions {x{t,to)}to parametrized 
with to- They are all on the exact curve W{t) t — to up to O(e^), but only vahd locally 
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for t near to. So it is conceivable that the envelope E of {Cto}to which contacts with 
each local solution at t = to will give a global solution. Thus the envelope function (t) 
coincides with W(t); 

x^{t) = x{t,t) = W{t). (3.13) 

Our task is actually to determine A{tQ) and ^(to) as functions of to so that the family of 
the local solutions has an envelope. According to the classical theory of envelopes given 
in the previous section, the above program can be achieved by imposing that the envelope 
equation 

= 0, (3.14) 

gives the solution to = t. From Eq.'s (3.11) and (3.14), we have 

dA ^ de e"^ , , 

where we have used the fact that dA/dt is 0(e) and neglected the terms of O(e^). Solving 
the equations, we have 

A{to) = Ae-^*«/^ ^(to) = -^to + e, (3.16) 

o 

where A and 9 are constant numbers. Thus we get 

x^{t) = x{t,t) = Woit) = Aexp(-^t) sin((l - ^)t + 9), (3.17) 



up to O(e^). Noting that y'l - 6^/4 = 1 - + 0{e^), one finds that the resultant 
envelope function ^^(t) = Wo{t) is an approximate but global solution to Eq.(3.1); see 
Eq. (3.2). 



4 Nonlinear equations 



In this section, we treat a couple of examples of systems of ODE's with nonlinearity to 
show how the RG method works[^]. The examples are the Lotka-Volterra|T^ and the 
Lorenz[|l4| equation. We shall derive the time dependence of the solution to the Lotka- 
Volterra equation explicitly; a phase equation will be derived by our method. The Lorenz 
equation is an example with three degrees of freedom, which shows a bifurcation. We 
shall give the center manifolds of this equation around the first bifurcation point. 
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4.1 Lotka-Volterra equation 



The Lotka-Volterra equation reads JIB 



X = ax — exy, y = —by + e'xy, (4.1) 

where the constants a, b, e and e' are assumed to be positive. It is well known that the 
equation has the conserved quantity, i.e., 

61n + a In \y\ — {ex + ey) = const.. (4.2) 



The fixed points are given hj {x = 0,y = 0) and (x = b/e',y = a/e). Shifting and 
scaling the variables by 

x={b + eO/e', y = a/e + r), (4.3) 

we get the reduced equation given by the system 

d 



where 



The eigen value equation 



has the solution 



(^-^o)u^-.e.(_\), (4.4) 



^o=f° ^. (4.5) 



7] J' ^ \a 

LqU = AoU (4.6) 



Let us try to extract the global behavior of the solution around the fixed point. Our 
strategy is the following: We suppose that we are on the exact solution at t = to where 
to is arbitrary; we denote the initial value by W(to)- We also suppose that we can apply 
perturbation theory for the solution at least in the small neighborhood of t = to- We 
expand the variable in a Taylor series of e; 

u = uo + eui + e^U2 H , (4.8) 

with Uj = An essential point of our method is to expand the initial value, too; 

W(to) = Wo(to) + eWi(to) + e=^W2(to) + ■ ■ ■ . (4.9) 

Our central task is to extract the initial value as a function of to so that the resulting 
local solutions starting from the different initial points at t = to and , say, t = to + At 
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are continued smoothly. This condition is found to be nothing but the one that the local 
solutions have an envelope. In actual calculations, it is also important to use the fact 
that the functional form of the initial values can be reduced from the general solution 
of the differential equations in a perturbative way; in this procedure, only independent 
functions modulo to secular terms are retained. 



The lowest term satisfies the equation 

d 



dt 



- Lq Uo = 0, 



which yields the solution 

Uo{t;to) ^ A{to)e''^^U + c.c. 
Notice that ^ is a complex number, so one may parametrize it as 

A{to) ^ A{to)/2i ■ eMi^ito))- 
The solution implies that the initial condition is given by 

Uo(to;to) =.4(to)e'"*°U + c.c.. 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



It means that in the lowest approximation the solution is parametrized by a complex 
function A(to) or a pair of real functions, A(to) and 0(to) = <^^o + ^(^o)- With a small 
perturbation, we expect that A and 9 will move slowly. 



Noting that 



-1 



aU + c.c, with a = (1 — ib/u!)/2, one finds that the first order 



term satisfies the equation 
/d 



iA^e'^^^aV + c.c.) + c.c. 



the solution to which may be given by 



A'iaV + — U*)e2^'^* + c.c. 
3 



(4.14) 



(4.15) 



Thus the initial value Wi(to) in this order is given Ui(to,^o)- 
Similarly, the second order solution may be given by 



U2 = 



h — iuj . 



+ 



362 



a 



\AYA\a{t-to + i—)\5 



2uj' 



2iuj 



-U* e 



+C.C. 



(4.16) 



Here we have a secular term proportional to the unperturbed solution. Since we want 
to make the lowest initial value as close as the exact one, we demand that as many as 
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possible terms in the higher order vanish at t = to- Thus adding unperturbative solutions, 
we make the secular term (of the upper component) vanishes at t = to-0 Thus neglecting 
higher order terms, we have u(t, to) = Uq + eui + e^U2, and W(i(:o) = u(to, ^o)- We impose 
that the solutions u(t, to) and u{t,tQ + At) give the same value at t. By taking a limit 
At — >• 0, we have the envelope equation; 

du 

with to = t. This gives the equation for A(t) as 



In terms of A{t) and 0{t), we have 



A(t) = const., 9{t) = -^-—{1 + —)u;t + 9o, (4.19) 

z4 uj 



with 6q being a constant. Owing to the prefactor i in r.h.s. of Eq. (4.18), the absolute 
value of the amplitude A becomes independent of t, while the phase 6 has a t-dependence. 
The envelope function is given by 

Ui.(t)=(^^[^^])=u(t,t)=W(t). (4.20) 

In terms of the components, one has 

C = Asine(t) - e— (sin2e(t) + ^cos2e(t)) 

6uj b 

€^A^3uj^-b\ . , Aub 

- -5[(^-^^)-^®W-T26^^-®^^) ^^-21) 



where 



2 ^2 

e(t)=cut + ^o, cD = {l-^(l + — )V. (4.22) 

24 cu^ 

One sees that the angular frequency is shifted. 

We can identify U£;(t) = {CE{t),VEit)) = W(t) as an approximate solution to Eq.(4.4) 
by construction. We see that u^(t) is an approximate but uniformly valid solution to the 

Although we can make other terms also vanish at i = to|ll]i the resulting dynamics of A and 6 
become more complicated than the present choice. In the theory of reduction of dynamics, one usually 
prefers simpler dynamics Pq| . 
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equation up to O(e^). We remark that the resultant trajectory is closed in conformity 
with the conservation law Eq. (4.2). 



"Explicit solutions" of two- pieces of Lotka-Volterra equation were considered by Frame 
r?! ; however, his main concern was on extracting the period of the solutions in an average 



method, we are not aware of any other work than ours which gives an explicit form of 
the solution as given here. 



4.2 The Lorenz model 



The Lorenz model [14] for the thermal convection is given by 



C = iv-K- (4.23) 

The steady states are give by 

(A) (e,r/,C) = (0,0,0), (B) (e, r/, C) = {±^h{r - 1), ±^h{r - 1), r - 1). (4.24) 



The linear stability analysis||T8[ shows that the origin is stable for < r < 1 but 
unstable for r > 1, while the latter steady states (B) are stable for 1 < r < a{a + h + 
3)/(cr — 6 — 1) = Tc but unstable for r > r^. In this report, we examine the non-linear 
stability around the origin for r ~ 1 and extract a center manifold for the dynamics in 
this region. 

We put 

r = n and ji = ■, X = sgn//. (4-25) 
We expand the quantities as Taylor series of e: 

u = 1^77 j = eui + e^ua + e^Ug + • • • , (4.26) 

where = *(^j, rji, d) {i = 1, 2, 3, . . .). We also expand the initial value at t = to; 

u(to, to) = W(to) = eWi + e^Wa + e^Wg + ■ ■ ■ . (4.27) 

The first order equation reads 

|-i„)u,=0, (4.28) 
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where 

/-a a \ 

Lo=\ I -1 , (4.29) 

V 0-6/ 

the eigenvalues of which are found to be 

Ai = 0, A2 = -cr - 1, A3 = -b. (4.30) 
The respective eigenvectors are 

f \ /o\ 

(4.31) 













■1 








) 













We are interested in the asymptotic state as f — > 00. In this asymptotic region, one 
may take only the neutrally stable solution 

Ui{t;to) = A{to)Vi, (4.32) 

because the other terms proportional to U2,3 will decay out at a sufficiently large time. 
Here we have made it explicit that the solution may depend on the initial time to- 
Eq.(4.32) implies that we have taken the initial condition that 

W(to)^Wi(to)=^(to)Ui. (4.33) 

In terms of the components, 

Ut) = A{to), Vi{t) = A{to), Ci{t) = 0. (4.34) 

In another word, the motion of Ui (t) is confined or reduced to the one-dimcnsional mani- 
fold '(A, A, 0), although A is a constant in this approximation. One expects that the small 
perturbation with e 7^ will give rise to a slow motion of A as well as a modification of 
the slow manifold. 

The second order equation reads 

/ 



Ui 1 = A^Us, (4.35) 



which may yield 



U2(t) = ^Us, (4.36) 



or in terms of the components 



6 = ?72 = 0, C2 = y. (4.37) 



11 



Here we have retained only functions independent of the ones appearing in the lowest 
approximation. Then the third order equation is given by 

^ - Lo)u3 = [ -x6 - ^Ci - eiC2 I = Y^^XA \A'){a\J, - U2), (4.38) 



6?7i + Ci?72 



which may yield 



U3 ^-^{xA- \A'){<j{t - to + T^)Ui - T^U^}. (4.39) 
l + cr 1 + (T 1 + cr 

Here we have again retained functions which have not appeared before except for the term 
with which the secular terms vanishes at t = to- Of course, one may have other choices 
for the independent functions, but it is found that the present "minimal" choice gives the 
simplest dynamics for the amplitude A. 

Thus collecting all the terms, one has 

n{t-to) = eA{to)\Ji + -A{tof\]^ 

+ ^ixA{to) - \A{tof ){o{t - to + TT^)Ui - T^U^}, (4.40) 

up to 0{e^). Accordingly, the initial value reads 

e2 

W(to) = u(to,to) = eA(to)Ui + -A(to)'U3 



,(x^(^o) - \A{tof){a\J, - U2}. (4.41) 



Demanding that the solutions at different initial times are continued smoothly, we 
have the RG/E equation, which reads 



dw. 
ato 



to=t 



^ e\xA-\A')\Ju 



l + (7 



(4.42) 



up to 0{e^). Noting that one may self-consistently assume that dA/dt — 0{e^), we have 
the amplitude equation 

dA 



,2 ^ I^.AI^\ 1 Al^\5 



(4.43) 



With this A{t), the envelope function is given by the initial value 

UE{t) = u{t;to = t), 
= W(t), 



2 3 

= eA{t)Vi + ^-A{tfVs+ ^ 



(l + <^) 
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^{xA{t)--A{tr){aV,-V,), (4.44) 



or 

Ut) = eA{t), r^j,{t)=eA{t) + -^{xA{t)-U{tf), (^(t) = '^A{ty . (4.45) 

1 + 0" 

We see that the initial value is obtained as the envelope of the local solutions and becomes 
a global solution to the Lorenz model. 

A remark is in order here; Eq.(4.45) shows that the slow manifold which may be 
identified with a center manifold |18[ is given by 

' ^ 1 + 0-^ b{l + a) b ^ ^ 

Notice that the invariant manifold is modified with the perturbation and also the slow 
dynamics Eq.(4.43) on the manifold is obtained by our method. One thus sees that 
the RG method is a powerful tool to extract center manifolds in a concrete form. It 
is worth mentioning that since the RG method utilizes neutrally stable solutions as the 
unperturbed ones, it is rather natural that the RG method can extract center manifolds 
when exist 0. 



5 The basis of the RG method for systems 

In this section, we give an account of our method in a general setting for ordinary equations 
as a summary. 

Let X = *(Xi, X2, ■ ■ ■ , Xn) and F(X, t; e) = *(Fi(X, t; e), F2(X, t; e), ■ ■ ■ , F„(X, t; e)), 
and X satisfy the equation 

^ = F(X,t;6), (5.1) 

with the initial condition X(to) = W, where to is arbitrary. We remark that the initial 
value W may be dependent on Iq, i.e., W = W(to)- We also write the solution X as 
X(t; to) so that the initial-time dependence is explicit. 

Let us try to have the perturbation solution of Eq.(5.1) around t = to by expanding 

X(t; to) = Xo(t; to) + eXi(t; to) + e'X2(t; to) ■ ■ ■ . (5.2) 

We also expand the initial value as 

X(to, to) = W(to) = Wo(to) + eWi(to)) + e'W2(to) + ■ ■ • . (5.3) 

In fact, Xj(to,to) = Wj(to). We suppose that an approximate solution X = X(t;to) to 
the equation up to O(e^) is obtained. It implies that the initial value W(to) at t = to 
coincides with an exact solution up to O(e^). We also have 

rfX(t;to^,W(to)) ^ p^^^^. ^^^^^^^ ^. ^ ^^^^^ 
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One may say that now we have a family of the orbits given by the functions X(t; to, W()f:o)) 
with to parameterizing the orbits. We see that the envelope E of the family of the orbits 
which contacts with each curve at t = to will give an approximate but global solution of 
the equation. Thus the envelope function is nothing but the initial value as a function of 
the initial time 

X^(t) =X(t;t,W(t)) = W(t). (5.5) 
The construction of E is performed as follows: We impose that the RG/E equation 



dX 

dto 







(5.6) 



gives the solution to = t, from which the dynamics of the initial value W(t) is obtained. 
Eq.(5.6) may give equations as many as n which are independent of each other. In the 
applications to describe asymptotic behavior of solutions, the equation is usually reduced 
to a low-dimensional equation. 

In accord with the above relation, one can easily show that X£;(t) = W(t) satisfies 
the original equation uniformly up to O(e^). In fact, Vto, one has 



dXf 



dt 



rfX(t;to,W(to)) 



t=to 



dt 

t/X(t;to,W(to)) 



rfX(t;to,W(to)) 



t=to 



dto 



t=to 



dt ^^^^ 
F(X£;(to),to;e) + oV) 



(5.7) 



where Eq.(5.6) has been used in the last equality. 



6 Concluding remarks 



We have described the perturbative RG method for global and asymptotic analysis. We 
have emphasized its relevance to the classical theory of envelopes and that the method 
concerns with the initial or boundary values of differential equations. It should be re- 
marked that the notion of envelopes is also useful for improvement of perturbation series 
in quantum filed theory^] and for obtaining asymptotic behavior of wave functions in 
quantum mechanics pO[] . 
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